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Abstract 

| ! Self-Organized Criticality (SOC) phenomena could have a signif- 

J> ' icant effect on the dynamics of ecosystems. The Bak-Sneppen (BS) 

■ model is a simple and robust model of biological evolution that exhibits 
, punctuated equilibrium behavior. Here we will introduce random ver- 
| sion of BS model. Also we generalize the single objective BS model 

to a multiobjective one. 

^ ■ Keyword: Self-organized criticality, evolution and extinction, BS model, 

multiobjective optimization. 

7? ■ 

■ 1 Introduction 

In recent years, there has been increasing interest of the possibility that 
evolution of species in an ecosystem may be a self organized critical phe- 
nomena. There are interactions among species in that ecosystem. The most 
common such interactions are predation, competition for resources and mu- 
tualism. As a result of these interactions the evolutionary adaptation of one 
species must affect its nearest neighbors. These interactions can give rise 
to large evolutionary disturbances, termed revolutionary avalanches. Most 
of these evolution models, like BS model [1], considered only one fitness for 
each species. Bak and Sneppen proposed a self organized model to explain 
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the punctuated equilibrium of biological evolution. They considered a 1- di- 
mensional model with periodic boundary conditions, topologically a circle. 
Assign a fitness < / (i) < 1 to each site i, i = 1,2,..., N, where N, is the 
number of species in the ecosystem. At each time step look for the site with 
lowest fitness j then replace its fitness together with the fitnesses of its near- 
est neighbors, j ± 1, by new ones which are uniformly distributed random 
variables. 

/ (j) = random value 
f (j + 1) = random value 
f (j — 1) = random value 

After running the system for sufficiently long time most of the fitness are 
above certain threshold .667. Also, the distribution of the distance between 
subsequent mutations and the avalanche sizes exhibit power laws. 

Several modification can be done to the BS model. The first possible 
modification is to use extremal dynamics [2, 3] which depends on the following 
idea: In real biological systems not only the lowest one who is updated but 
some of the low fitness species. This number that changes is not fixed but 
random. So, we will study this random version of BS model in section 2. 
Also, in biology almost every optimization problem is multiobjective (MOB) 
e.g. objective of foraging and of minimizing predation risk. In section 3, we 
will apply the concept of MOB to BS model. 

2 Random BS Model 

Here, we will study the first modification that can be done to the BS model. 
Instead of finding exactly the site with lowest fitness, one may use the ex- 
tremal dynamics. In this case a uniformly distributed random number is 
picked and all the sites with fitness less than this number has its fitness up- 
dated. This dynamics has been used to explain the long term memory for 
the immune system [3]. It has been also used to solve some optimization 
problems e.g. spin glass, graph coloring and graph partitioning [4, 5]. We 
run a system consisting of N = 4096 species for different sufficiently long 
time (up to 2 x 10 7 ). We find that most of the fitnesses are above a certain 
threshold value 0.64, as shown in figure 2. In figure 1 we plotted the standard 
Bak-Sneppen model for reference. 
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3 Multiobjective Optimization Model 



In most evolution models e.g. BS model, only one fitness is considered i.e. 
single objective optimization. Almost every real life problem is multiobjective 
(MOB) one [6]. Therefore it is important to generalize the standard single 
goal oligopoly studies to multiobjective ones. Methods for MOB optimization 
are mostly intuitive. 

The first method is lexicographic method. In this method objectives are 
ordered according to their importance. Then the first objective is satisfied 
fully. The second one is satisfied as much as possible given that the first 
objective has already been satisfied and so on. A famous application is in 
university admittance where students with highest grades are allowed in any 
college they choose. The second best group are allowed only the remaining 
places and so on. This method is useful but in some cases it is not applicable. 

The second method is the method of weights [7]. Assume that it is 
required to minimize the objectives Z(i), i = 1,2,..., AT. The problem of 
maximization is obtained via replacing Z(i) by —Z(i). Define 

N 

z = y: w(i) z(i) 
i=i 

where 

N 

w(i) > and ^2w(i) = 1 

i=i 

Then the problem becomes to minimize Z. This method is easy to implement 
but it has several weaknesses. The first is that it may give a Pareto dominated 
solution. A solution Z (i), i = 1,2, ...,N is Pareto dominated if there is 
another solution Z(i), i = 1,2, N such that Z{%) < Z (i) for all i with at 
least one k such that Z(k) < Z (k). The second difficulty of this method is 
that it is difficult to apply for large N. 

The third method is to minimize only one objective while setting the 
other objectives as constraints e.g. minimize Z{1) subject to Z(i) < a(i), 
% = 2,3, ...,N where a(i) are parameters to be updated. The problem with 
this method is the choice of the thresholds a(i). In the case of equality i.e. 
Z{%) = a(i) this method is guaranteed to give a Pareto optimal solution. 

The fourth method using fuzzy logic is to study each objective individ- 
ually and find its maximum and minimum say Z max (z), Z min (i) respectively. 
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Then determine a membership 

Z(i)-Z max (i) 
mW = Z (?) - Z (j) 

Thus < m{i) < 1. Then apply max (min (m(i), i — 1,2, A/")). Again this 
method is guaranteed to give a Pareto optimal solution. This method is a 
bit difficult to apply for large number of objectives. 

The BS model can be generalized to the multiobjective. Assigning two 
fitnesses fi (i), f 2 (i), to each site instead of one. The updating rule is If 

x fi (i) + (1 — f 2 (i) < min fit 

where < x < 1, then update both fi (i) , f 2 (i) and fi (i ± 1), f 2 (i ± 1). 

In the updating rule we have used the simple and widely used method, 
weighting method in MOB. Multiobjective optimization is much more realis- 
tic than single objective ones. After running a system consisting of A^ = 4096 
species for different sufficiently long time (up to 2 x 10 7 ). The distribution 
of the distance between subsequent mutations are shown in figure 3. We find 
that most of the fitnesses are above certain threshold value 0.57, as shown in 
figure 4. The size of avalanches are shown in figure 5. 
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FIGURES 



FIG. 1. The results of BS model under our simulations in a system of size 
N = 4096 and t = 10 7 , iteration, (a) The distribution of distances 
D(x) between subsequent mutations x. (b) Distribution of fitness in 
the critical state (right curve) D(F) with the distribution of minimum 
fitness (left curve), (c) Distribution of avalanche sizes D(S) in critical 
state, (d) Mutation activity vs time measured as the total number of 
mutations (N = 64 and t = 4x 10 3 ). 

FIG. 2. The results of random BS model in a system of size N = 4096 and t = 
10 7 , iteration, (a) Distribution of fitness in the critical state D(F). (b) 
Distribution of minimum fitness in the critical state, (c) Distribution 
of avalanche sizes D(S) in critical state. 

FIG. 3. The distribution of distances D(x) between subsequent mutations in a 
system of size iV = 4096 and t = 10 7 , iteration with three different 
weights 0.3, 0.5, 0.9. 

FIG. 4. The distribution of fitness in the critical state (right curve) D(F) with 
the distribution of minimum fitness (left curve) in a system of size 
N = 4096 and t = 10 7 , iteration with three different weights 0.3, 0.5, 
0.9. 

FIG. 5. The distribution of avalanche sizes D(S) in critical state in a system of 
size N = 4096 and t = 10 7 , iteration with three different weights 0.3, 
0.5, 0.9. 
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Figure 1 
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Figure 5 



